Conformational Instability of Rodlike Polyelectrolytes due to Counterion Fluctuations 
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The effective elasticity of highly charged stiff polyelectrolytes is studied in the presence of coun- 
terions, with and without added salt. The rigid polymer conformations may become unstable due 
to an effective attraction induced by counterion density fluctuations. Instabilities at the longest, 
or intermediate length scales may signal collapse to globule, or necklace states, respectively. In the 
presence of added-salt, a generalized electrostatic persistence length is obtained, which has a non- 
trivial dependence on the Debye screening length. It is also found that the onset of conformational 
instability is a re-entrant phenomenon as a function of polyelectrolyte length for the unscreened 
case, and the Debye length or salt concentration for the screened case. This may be relevant in 
understanding the experimentally observed re-entrant condensation of DNA. 



I. INTRODUCTION AND SUMMARY 

A polyelectrolyte (PE) is an ionic polymer which when 
dissolved in polar solvents dissociates into a long poly- 
mer chain (macroion), and small mobile counterions. Be- 
cause of the electrostatic repulsion of the uncompensated 
charges on the polymer, the chain is stretched out to rod- 
like conformations. Naively, one may expect that higher 
charges on the polyelectrolyte lead to more rodlike con- 
figurations. This is not the the tendency to- 
wards extended shapes is opposed by stronger attraction 
to the counterions. The latter may condense on a highly 
charged polyelectrolyte, giving it a much lower apparent 
charge and resulting in a lower stiS^ness. Nonethe- 
less, in a mean-field (Poisson-Boltzmann) treatment of 
the counterions, the conformation of the polyelectrolyte 
is stretched [|,|. 

On the other hand, it is now well known from ex- 
periments [^,||, simulations 0, and analytical theories 
|p[-p^ , that highly charged flexible polyelectrolytes (such 
as polystyrene sulphonate) can collapse in the presence 
of multivalent counterions to highly compact states. Stiff 
polyelectrolytes (such as DNA) could also collapse [11 1, 
and rodlike polyelectrolytes may form bundles [p^l4|, 
under similar conditions. To account for these effects, 
an attractive interaction capable of competing with the 
residual coulomb repulsion (of the polyelectrolyte back- 
bones with compensated charge densities) is needed. 
Note that the situation is more subtle for stiff polyelec- 
trolytes due to their intrinsic rigidity ll^Jl^ , [l8| -po| . 

To understand the origin of the attractive interaction, 
consider two like-charged substrates on which counteri- 
ons are condensed. As the two objects approach each 
other, the counterions may rearrange their positions. 
Any correlated separation of charges now leads to an at- 
traction whose range is of the order of the "correlation 



hole" of the counterions on each substrate. Several dif- 
ferent mechanisms may lead to correlated charge sepa- 
ration: (i) At low temperatures, the counterions form 
a Wigner crystal on the charged substrate, which leads 
to an attractive interaction with a range set by the lat- 
tice spacing pl|-p5|. (ii) Specific binding of counterions 
to the substrate is another mechanism for charge separa- 
tion H^] . In this case the periodicity, and hence the range 
of the attraction, is dictated by the underlying structure 
of the substrate, such as the double helical structure of 
DNA along which the charged phosphate groups are lo- 
cated, (iii) Thermal fiuctuations can also induce instan- 
taneous charge separations which inter-correlate on the 
two objects, leading to an attraction similar to the van 
der Waals interaction. Since the dominant fluctuations 
have wavelengths of the order of the separation of the 
objects, the range of the attraction is set by the dis- 
tance between them; i.e. the interaction is long-ranged 
§,|lj2|j2|^. (iv) Finally, it has been proposed that 
macroions can be overcharged by the conden sing coun- 
terions due to the gain in correlation energy [ pO| . Con- 
sequently, when two such macroions are brought close 
to each other in a neutralizing solution, the distribu- 
tion of the counterions can be asymmetric leading, at 
each instant, to one overcharged and one undercharged 
macroion ]3l| , |3^ . This mechanism will lead to instanta- 
neous ordinary long-ranged Coulomb attraction between 
the two decorated macroions, and is reminiscent of cova- 
lent bonding in atomic systems in that the macroions are 
effectively sharing a fraction of their condensed counte- 
rions I p^ . 

As thermal fluctuations clearly set the energy scale for 
the relevance of these different attractive interactions, 
one may naturally ask the following question ||3^,Q: Is 
room temperature 'high' or 'low' for DNA or actin con- 
densation, or bundle formation? In regards to this ques- 
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tion one should note that although an actual Wigner 
crystal may appear far-fetched on the surface of DNA 
or actin at room temperature, remnants of such a struc- 
ture may survive for multivalent counterions. 

Since it is just the local charge separation that is 
needed for the 'zero temperature' mechanism, and not 
the long-range order of the global structure, a liquid 
phase should still lead to attraction provided that some 
short-range order persists. Consequently, it seems that 
both types of interactions are in effect at room tempera- 
ture and the dominant mechanism should be determined 
by comparing the range of the interaction and the dis- 
tance between macroions. While the short-ranged inter- 
actions seem to play the dominant role in bundle forma- 
tion |2l|- [2^ , ^^ , there seems to be a consensus that they 
are not capable of actually collapsing a single DNA or 
actin chain , and thus the long-ranged interaction 

caused by thermal fluctuations may be responsible for 
the observed condensation |pl] , |l8| . 

In this paper, we focus on the long wavelength behavior 
of polyelectrolytes. Our approach is complementary to 
that of the "Wigner liquid" theories [|l6|,0 , whose local 
formulation corresponds to examination of short wave- 
length behavior. They find that the softening contribu- 
tion to the persistence length due to the Wigner liquid 
correlations is not large enough to fully compensate the 
mechanical stiffness and render a null persistence length, 
which would trigger chain collapse. We also find that 
the softening contribution to the rigidity due to charge 
fluctuations at the smallest length-scales is never strong 
enough to completely negate the bare bending rigidity, 
and obtain collapse at these scales, but rather that the 
'collapse' occurs at higher length-scales due to longer 
wavelength charge fluctuations. Additional correlations 
at short length-scales as might be obtained with a Wigner 
liquid should not affect the large-scale behavior which is 
the domain of our theory. 

To this end, we employ path integral methods [ p9|j35| 
to study the energy cost of deforming a stiff and highly 
charged polyelectrolyte in the presence of thermally fluc- 
tuating counterions [l^ . In particular, consider a chain 
of length L, with a microscopic persistence length £p, and 
average separation a between charges on its backbone, in 
a neutralizing solution of counterions of valence z. The 
polyelectrolyte is considered highly charged when a is less 
than the Bjerrum length £b = e^/^ksT, where e is the 
dielectric constant of the solvent. (For water at room 
temperature, £b — 7.1 A.) We calculate the effective free 
energy of a fluctuating polyelectrolyte as a perturbative 
expansion in the deformations around an average rodlike 
structure. The linear stability of this structure is con- 
trolled by a spectrum £{k), as a function of the defor- 
mation wave vector k. A negative value of £{k) signals 
an instability at the corresponding wavelength, leading 
to phase diagrams as in Figs. |^ and |^. In particular, 
both in the presence of salt and in salt-free conditions, 
(see Figs. || and ||), we find that counterion fluctuations 
cannot trigger collapse of a stiff polyelectrolyte, unless it 



has a microscopic persistence length less than a critical 
value of 

The index i refers to either added-salt (as) or salt-free 
(s/) conditions, with corresponding numerical constants 
Aas and Agj which are given later. For ip < ip, there 
is a finite domain of intermediate polyelectrolyte lengths 
L, for which a collapsed conformation is favored. 

For a low concentration of added salt, we find an effec- 
tive persistence length of 

^ ^ aia/iBf C2 

16z^(l-^)^.a)Mn^(i,)"2.1n(J,)' 

(2) 

where is the Debye screening length, related to the 
salt (number) density n via — AniBn-, and C2 is 
a numerical constant given below. The above expres- 
sion, which is a non-trivial generalization of the Odijk- 
Skolnick-Fixman (OSF) electrostatic persistence length 
p6| , is a sum of repulsive (-I-) and attractive (-) electro- 
static contributions. It can therefore be negative under 
certain conditions, which we take as an indication of a 
conformational instability of the rodlike polyelectrolyte 
(tendency to collapse) . This occurs in the regions of the 
phase diagram indicated in Fig. ^, and only for persis- 
tence lengths less than a critical value given by Eq. (|l|) 
above. A similar 'softening' contribution to the rigidity 
and an associated conformational instability, can also be 
obtained for charged membranes There are related 
instabilities caused by surface fluctuation-induced inter- 
actions between stiff polymers on membranes 

We also find that the onset of conformational insta- 
bility is a re-entrant phenomenon as a function of poly- 
electrolyte length for the unscreened case, and the Debye 
length or salt concentration for the screened case. This 
may be relevant in understanding the re-entrant conden- 
sation of DNA that has been observed experimentally 
p9[ , but note that alternative explanations for this phe- 
nomenon also exist in the literature 



II. THE MODEL 

We consider a polyelectrolyte oriented along the x-axis, 
and parameterize is transverse deviations from a straight 
line by a two component vector r{x). By assuming a sin- 
gle valued function r(x), we implicitly neglect overhangs 
and knots, as sketched in Fig. |l|. The embedding of the 
polyelectrolyte in space is thus described by the three 
component vector R{x) = (x,r(x)). For each conforma- 
tion of the polyelectrolyte, we would like to calculate a 
constrained partition function Z[R], by integrating over 
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all possible configurations of the counterions in the so- 
lution. We do this by using a phenomenological model, 
which generalizes previous work by Ha and Liu (and 
is referred to as the Ha-Liu model henceforth), within 
which the restricted partition function is written as 



Z[R] EE e = e X 



dx [q{x) 



X exp 
2 



1 



go 



Vq{x) 
2 



27 a 

dx dx' q{x)q 


{5qY 


a a 


R{x) - R{x') 





(3) 



In this equation, the energy cost for the deformations of 
the polyelectrolyte is characterized by an intrinsic bend- 
ing energy described by the Hamiltonian 



mm = f / dx{diRf 



(4) 



and a fluctuating electrostatic self energy. The fluctu- 
ating linear charge distribution on the PE is character- 
ized by a Gaussian distribution, with the mean (cor- 
responding to the average renormalized charge density 
and read off from the exact Fuoss-Katchalsky-Lifson so- 
lution for an infinitely long uniformly charged rod 41 1) 
given by go = oJzIb, and the variance given by {5qY = 
z{l — a/zls) It should be noted that we have in- 

cluded the effect of a low density of added salt in the 
above equation by adding a screening term The 
salt- free case corresponds to k = 0. 



III. EFFECTIVE ELASTICITY OF THE 
POLYELECTROLYTE 

We next use the methods of Refs. p9| , |35[ |, to calculate 
an effective Hamiltonian from the above Gaussian path 
integral, perturbatively in the deformation field r(a;). We 
can rewrite the expression in Eq. (0) above as 



Z[R] = e-''^''[«l-i^ X j Vq{x) 

X exp(-i I ——q{x)M{x,x')q{x') 
I 2 J a a 

^JSqj^J-'^^\ 



where 



M{x,x')^j^S{x-x') + '''' 



[dqf 



R{x) - R{x'] 



(5) 



(6) 



Functional integration over q{x) then yields 



/3H°ff = pnP[R] + TT^T^T + ^lndet{M(a;,a;')} 
Za(dq)^ 2 



1 On f dx dx' , , 1 , 

/ M-\x,x'). 

(0(7)4 J a a 



(7) 
(8) 



2 {5qY 
Now we can write 

M{x,x') = Mq{x - x') + 5M{x,x'), 
or alternatively in Fourier space 

M(fc, k') = 2TrS{k + k')Mo{k) + 5M{k, k'), (9) 
in which 



Mo{k) 



Ib ln[(fc^ + K2)a2i 



{5qf 



(10) 



and expand Eq. (Q) in powers of 5M. To the leading 
order, this yields 



= m^ + /3^o + ^ tr{M-^6M) 
1 ql 5M{k ^0,k' ^ 0) 



-t- 



2 {5q)* [Mo{k = 0)]2 



(11) 



In the above equation, the trace term is an attractive 
fluctuation-induced effective free energy, and the term 
proportional to gg is a repulsive electrostatic (self-) free 
energy. The constant term pj-Q is independent of the 
deformations, and thus plays no role in the effective elas- 
ticity of the chain. It will be neglected henceforth. 

Using the deflnition of M [Eq. (||)], and the expansion 
about the rodlike conflguration, we find to the leading 
order 



ln[(fc2 + K^)a^] 



+^[ip + fc) + « ] Hip + kya' + i^ay] 

+ ^-f[{p-ky + K^] Hip ~kfa' + [naf] | . (12) 

We can then use the expressions for Afo(fc) [Eq. (px| )] 
and 5M{k, -k) [Eq. (^] and put them into Eq. (|11|) to 
find the effective Hamiltonian. We find 
1 r i-ii. 

m^'^'irix)] - 2 y ^ ^(^) 1^(^)1' + (13) 
where 

i^Jxi^ — - 

Qo^B { (k^ + K^) ln[(/c2 + K2)a2] _ ^2 in[(^a)2] 1 



2a2 



[l + 2(if)(5g)2 1n(J,) 



^ 2 /■ dp I (p^ + K^)ln[(p^ + ^V] 
2a^ J 2tt\i- (^) (,5g)2 ln[(p2 + ^2)02] 



[jk + pf + K^] ln[(fc + pfa^ + jnaf 
1- (^) ((5g)2ln[(p2 + K2)Q2] 



+ tP. (14) 
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The above expression for the elastic kernel consists of 
four different contributions: 

(i) The intrinsic rigidity of the polyelectrolyte, which 
has a regular dependence corresponding to curvature 
elasticity. 

(ii) A repulsive contribution corresponding to electro- 
static stiffening of the polyelectrolyte. This is a gener- 
alization of the wave vector dependent rigidity first in- 
troduced by Barrat and Joanny and has a non- 
analytical k dependence for k = 0. (In fact, one can 
simply recover their formula by setting {Sq)"^ — 0.) 

(iii) A fluctuation-induced attractive term that com- 
petes with the previous two contributions and tends to 
soften the polyelectrolyte. 

(iv) A term corresponding to line tension renormaliza- 
tion with a k^ dependence. In general, one can show that 
such a term should be absent since the rotational sym- 
metry of the original Hamiltonian for a polyelectrolyte 
requires £{k)/k'^ ^ in the hmit fc ^ |Q. Since the 
other terms in Eq. ( p^ ) above have a non- vanishing limit 
as fc — >• 0, one can simply tune r in such a way that the 
overall line tension vanishes in that limit |Q . 

We now examine two limiting cases of the above result 
in more details. 



A. Salt-free solution 

The salt-free case is relevant to situations in which the 
screening length is much larger than the length of the 
polyelectrolyte, i.e., for ^ L. Also in this limit the 
deformation energy is a non-analytic function of fc, and 
taking the k — limit of Eq. ([l^ ) requires some care. 



{i)^B fc2ln[(A)Vl 



2[l + 2W(^)ln(^)]" 



-B^\k\ 



(15) 



with fco <C IT / L being a cutoff which in the absence of 
screening is required for a polyelectrolyte ^] . Note that 
i9fe2f (fc)|fc_>o = to ensure the rotational symmetry, and 
that Eq. (|l^) is not just the first few terms in a low-/c 
expansion. The constant Bi is 



Bi 



L/a^^(l+x2) 



In 



2 

X 


- 2a; In 


x + l 




x-1 





An 



(a/iB) 



ln(^ 



(16) 



The integral Bi very slowly depends on {5qf'{tB / a) for 
typical values, and can be best approximated as Bi ~ 
ci/ln^ (f) with ci ~ 0.101. 

The spectrum of the deformation modes in Eq. (|1E 
consists of an electrostatic repulsion term (that yields 
the familiar spectrum of a flexible polyelectrolyte in the 
limit {5qY = @|7))- ^"^^ fluctuation-induced attrac- 
tion reduces the energy cost, and could even lead to an 



instability at shorter length scales |Q . The rigidity term 
Ipk'^ is the largest power of k included in Eq. (|l5| ) and 
ensures stability at the shortest scales. A negative value 
of £{k) for any k indicates a linear instability; the on- 
set of which can be tracked by finding the point when 
when the minimum of the spectrum, determined from 
d£{k)/dk = 0, hits the line £{k) — 0. This corresponds 
to the criterion £p = i^, where the critical persistence 
length £p is given by Eq. (|^) above, with 



(17) 



In^ 



2£pfco ln^(L/a) 



For £p < there is domain of unstable modes for 
A:_ < k < k+, where 



Cl 



2ep\n{L/a) 



l±Jl 



(18) 



The spectrum of modes given by Eq. (|T^) is plotted in 
Fig. IJ. To determine the effect of the unstable modes 
on the conformation of the polyelectrolyte, we should 
compare the wave vectors with tt/L. There are three 
possibilities: (i) For fc+ < tt/L, the unstable modes can- 
not be accessed and the polyelectrolyte has an extended 
structure, (ii) For k- < tt/L < fc_|_, long wavelength 
modes have negative energy. Their unstable growth is 
likely to lead to the collapse of the whole chain, (iii) For 
tt/L < k-, the longest wavelength t^/L) are stable, 
but there is a range of unstable wavelengths. While lin- 
ear stability cannot determine the eventual state of the 
polyelectrolyte, the presence of a finite length scale could 
well be the precursor of a necklace structure. 

The corresponding phase diagram is depicted in Fig. ||. 
A suggested boundary between the collapsed and other 
phases is obtained from k± — tt/L, as 



In 



Cl 



Lkr 



(19) 



and has a nearly parabolic shape, leading to reentrant 
extended states. 



B. Added-salt 

The addition of salt leads to a finite screening length 
and for <C L, the spectrum £{k) given by 
Eq. (|l^) is an analytic function of k with a well defined 
power expansion. The lowest order term in the expansion 
is proportional to k^, whose coefficient can be regarded 
as the 'effective rigidity' of the polyelectrolyte, and given 
by 



L„ — 



1 



El 
2k' 



(20) 
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in which 



Bo 



1 / na 



dx 



llx'^ + 2x^ - 1 



An 



In 



(1/ Ka 
TT 



' + 1 



(21) 



Again, the integral B2 has a very slow dependence on 
{Sq)^{£B/a) for typical values, and can be best approx- 
imated as B2 — C2/ln(-^) with C2 — 0.288, leading to 
the effective persistence length given in Eq. (^. It is in- 
teresting to note that a similar result has been predicted 
for stiff polyampholyt es |pO[ ]. 

Interestingly, Eq. ( po|) reproduces the OSF electro- 
static persistence length in the limit (Sq)^ = |^^, with 
a reduced charge density qq. Upon including the coun- 
terion fluctuations, (Sq)'^ / 0, there is a reduction of 
the repulsive term, as well as the appearance of an at- 
tractive term. As a result, the polyelectrolyte 'rigidity' 
(effective persistence length) can become very low and 
even have negative values which we take as indicating 
a conformational instability (collapse). Using this cri- 
terion (Lp = 0), we obtain the phase diagram shown in 
Fig. 4, which like the salt-free case, has a nearly parabolic 
shape, with a maximum that yields the critical persis- 
tence length given in Eq. (0) with Aas = 4 ^ 0.0829. 
Once more, its shape suggests that the instability is a re- 
entrant effect. Note that the critical persistence length 
for the case with added salt is larger than the salt-free 
case. 



IV. DISCUSSION 

The above results can be related to some experimen- 
tal findings concerning polyelectrolyte elasticity in the 
presence of multivalent counterions. In a direct single 
DNA manipulation experiment, Baumann et al. used 
a force-measuring optical tweezer to determine the elas- 
tic properties of A-bacteriophage DNA in the presence 
of multivalent counterions They observed that the 
apparent persistent length of the DNA (extracted from 
force-extension curves using the wormlike chain model) 
goes down to as low as 250-300 A, which is well below 
the fully saturated high salt value 450-500 A (denoted as 
(.p above) in the presence of multivalent ions ||5l|. This 
is in agreement with the generalised effective persistence 
length given in Eq. (^. There may be of course, other 
contributions to the effective persistence length coming 
from other types of charge correlations. 

Flexible polyelectrolytes have also been shown to col- 
lapse in the presence of divalent ions |^]. Interestingly, 
there is evidence that for sodium polyacrylate chains the 
collapsed states are not always compact spherical shapes 
at lower salt concentrations and may have cigar-like 
or pearl-necklace shapes. This is in agreement with the 
results we obtained in the salt-free regime. 

We now conclude with a qualitative summary of the 
nature of the results, and the range of their validity. From 



dimensional analysis, it is easy to show that unscreened 
Coulomb interactions make a contribution of (.B{k/ to 
the rigidity spectrum £{k). In the PB solution, due to 
charge condensation the strength of this term is reduced 
by a factor of q^ — iajzlBf' ■ If the charge density on 
the polyelectrolyte is allowed to fluctuate, it is further re- 
duced, and the Coulomb rigidity goes down by a factor of 
{8qf{lBlaf\^{Lla), with {6qf = z(l - qo). However, 
these reductions do not change the overall sign which still 
prefers a rodlike structure. An attractive (destabilizing) 
contribution is generated by fluctuation-induced interac- 
tions, which are typically independent of microscopic pa- 
rameters, and hence make a contribution of —k^ to £(k). 
In comparison to the leading Coulomb contribution, the 
latter corrections become important at short scales of or- 
der of (P' IIb- We thus may well question the applicabil- 
ity of continuum formulations to describe such a short- 
distance instability. In hindsight, the phase diagrams of 
Fig. 2 and Fig. 3 indicate that the prefactors involved in 
softening of the residual repulsion conspire to make the 
actual instability lengths quite large (~ z^^^/a^), and 
thus the continuum formulation should hold for a large 
portion of these phase diagrams js^ . 

Another potential concern is the choice of the charge 
variance ((5<?)^. While our approach ignores the finite size 
of the counterions and allows for a large number of them 
to be condensed in the vicinity of a single charged group 
on the polyelectrolyte, there are other models that have 
restricted this number to one |^,^ . Given the relatively 
large radius of DNA and actin as compared to the size 
of typical counterions, the restriction of this number to 
one seems to be artificial, as it leads erroneously to zero 
charge fluctuations at complete condensation. Nonethe- 
less, we believe that considering the finite size of the 
counterions in a realistic way will certainly be a worth- 
while improvement of our approach. Finally, the instabil- 
ity analysis performed here only provides us with infor- 
mation concerning the onset of a conformational change. 
The final structure of the collapsed chain is naturally be- 
yond this linear stability analysis [ pTp^ . 

It is known experimentally that different counterions 
with the same valence may behave differently as collaps- 
ing agents. The difference is usually attributed to mod- 
ifications in the microscopic structure of the polyelec- 
trolyte that take place upon binding of the counterions 
It is thus plausible that the microscopic features 
that distinguish between different counterions with the 
same electrostatic properties can be encoded in a sin- 
gle parameter, the microscopic persistence length, which 
normally depends only on the local microscopic structure 
of the polyelectrolyte backbone. 
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FIG. 3. Phase diagram for the added-salt case, for different 
values of the counterion valence z. 




FIG. 1. Schematic of the extended polyelectrolyte 
parametrized by R{x) — {x,r{x)). 



7 



FIG. 4. The energy spectrum £{k) for the salt-free regime 
showing the unstable modes. 
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